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Abstract 

The three species asymmetric ABC model was initially defined on 
a ring by Evans, Kafri, Koduvely, and Mukamel, and the weakly asym- 
metric version was later studied by Clincy, Derrida, and Evans. Here 
the latter model is studied on a one-dimensional lattice of N sites with 
closed (zero flux) boundaries. In this geometry the local particle con- 
serving dynamics satisfies detailed balance with respect to a canonical 
Gibbs measure with long range asymmetric pair interactions. This 
generalizes results for the ring case, where detailed balance holds, and 
in fact the steady state measure is known only for the case of equal 
densities of the different species: in the latter case the stationary states 
of the system on a ring and on an interval are the same. We prove that 
in the limit N ^ oo the scaled density profiles are given by (pieces of) 
the periodic trajectory of a particle moving in a quartic confining po- 
tential. We further prove uniqueness of the profiles, i.e., the existence 
of a single phase, in all regions of the parameter space (of average den- 
sities and temperature) except at low temperature with all densities 
equal; in this case a continuum of phases, differing by translation, co- 
exist. The results for the equal density case apply also to the system 
on the ring, and there extend results of Clincy et al. 
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1 Introduction 

One dimensional systems play an important role in statistical mechanics. In 
addition to their intrinsic interest as models of physical systems in confined 
quasi-linear geometries, they are in many cases exactly solvable and in fact 
are, with few exceptions, the only exactly solvable many body systems 111 
[21 [3l H] . These exact solutions provide insights into both equilibrium and 
non-equilibrium collective behavior in higher dimensions. 

An interesting connection between equilibrium and non-equilibrium 
phase transitions in one dimension is provided by the ABC model intro- 
duced by Evans et al. [5j (a model with similar behavior was discussed in 
[6]): a one dimensional system consisting of three species of particles, labeled 
A,B,C, on a ring containing N lattice sites (one may equivalently regard it 
as a model with two species and empty sites, or as a "clock" model with three 
states). We will here typically let a = A, B , or C denote a particle type, and 
make the convention that a + 1, a + 2, ... denote the particle types which 
are successors to a in the cyclic order ABC. The system evolves by particle 
conserving nearest neighbor exchanges with asymmetric rates: 07 — > 7a 
(clockwise) with rate qa,^, where qa,a+i < Qa+i,a', the total numbers Na 
of particles of each species are conserved and satisfy J2a = Evans 
et al. argued that this system will, in the limit N ^ 00 with N^/N — > 
and with Tq, > for all q, segregate into pure A, B, and C regions, with 
rotationally invariant distribution of the phase boundaries. For the case 
QAB = QBC = QCA = q < 1 and qba = QCB = QAC = 1, the only case we 
will consider here, they showed further that when A^^ = A'^^ = = 
the dynamics satisfies detailed balance with respect to the Gibbs measure 
of a certain Hamiltonian having a long range pair interaction, so that the 
stationary state is an equilibrium state. 

In a later development Clincy et al. [7J considered a weakly asymmet- 
ric version of this model in which q = e~^^^; the stationary state for the 
equal density case = N/3 then becomes a Gibbs measure of the form 
exp{— /3£^7v}i with energy Ejy equal, up to a constant, to the Hamiltonian of 
[5] divided by A^. The parameter P thus plays the role of an inverse temper- 
ature: T = (The energy is written as En in [7j, but here we adopt the 
convention that quantities with and without a circumflex refer respectively 
to the ring geometry and to the interval geometry introduced below.) Clincy 
et al. obtained the Euler-Lagrange equations for the minimizers of the free 
energy functional, which here is equivalent to the large deviation functional 
(LDF), of the rescaled particle densities Pa{x), in the limit A^ — > 00 with 
Na/N — > 1/3 and i/N — > x, and showed that the uniform density profiles 
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Paix) = Ta = 1/3 are always a solution of these equations, but that for 
temperatures T below the critical temperature Tc = = {2it^/2>)~^ there 
is also a nonuniform solution of the linearized equations which has a lower 
free energy than the uniform state. They interpreted this as a second order 
transition at Tc from the uniform to the nonuniform state, and confirmed 
by numerical simulations that for T > Tc the system is in a single phase 
with essentially no correlations (they are 0{1/N)) between the locations of 
particles of different species, while for T < Tc there is segregation of the 
different species. They further argued that this transition will persist (pos- 
sibly becoming first order) for unequal densities. In this case the stationary 
state is no longer an equilibrium state and is in fact unknown; despite this, 
they were able to obtain the LDF to order for all densities (see also [8j). 

In the present work we consider the weakly asymmetric ABC system 
on a one-dimensional lattice of sites i = 1,...,N with zero flux boundary 
conditions: the dynamics are the same as above, except that a particle at 
site i = 1 (respectively i = N) can only jump to the right (respectively 
left). We shall refer to this geometry as an interval. In contrast to the 
situation on the ring, there is for this system always, whatever the values 
of Na, Nb, and Nq, an energy function En (see (j2.ip below) such that the 
dynamics satisfies detailed balance with respect to an equilibrium measure 
Up = Z-^ exp{-/?£'Ar}. 

When A^^ = Nb = Nc = N/3 the energy E]\[ for the system on the 
interval agrees with the energy En' on the ring [3, in the sense that if 
we (mentally) connect site N to site 1 clockwise, and thus identify each 
configuration in the interval with a corresponding configuration on the 

ring, then En{C) and i?7v(C) agree up to a constant which depends only on 
the Nee- Thus the probabilities of the configurations C and C (at the same 
/3) are the same, and the invariance under rotations of En implies a rather 
surprising "rotation" invariance of E^ and of the Gibbs measure on the 
interval. It follows also that all the results obtained in the present paper for 
the interval model in the case of equal densities give corresponding results 
on the ring. When the Nq. are not all equal one might of course use E^ 
similarly to define an energy of a ring configuration, but this energy would 
then depend on which site is chosen for the origin. The resulting Gibbs 
measure on the ring would be neither rotationally invariant nor invariant 
under the time evolution of the ABC dynamics. 

Using the Gibbsian nature of the invariant measure and following [7] we 
can obtain the free energy functional ^{pA, Pb, Pc) of the density profiles 
Paix) in the scaling limit A^ ^ oo, i/N x. The parameters of the model 
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in this limit are the inverse temperature (3 and the mean densities = 
Jq Pa dx. Our main goal is to determine, for given values of these parameters, 
whether has a unique minimizing profile or whether, conversely, there is 
more than one minimizer and thus coexistence of phases, and in either case 
the form of the minimizing profile(s). We are able to establish uniqueness at 
high temperatures by a direct study of the functional J-'. We further show 
that minimizers must satisfy the Euler-Lagrange equations for J^, and by 
a study of the solutions of these equations we establish the nonuniqueness 
of minimizers when = = rc = 1/3 and T < T^, the nonuniqueness 
corresponds to the lack of rotational invariance for the minimizers on the 
ring. Conversely, we establish uniqueness of the minimizers (on the interval) 
whenever the densities are not all equal. In the process we also obtain the 
form of the minimizing density profiles as elliptic functions corresponding to 
pieces of the periodic trajectory of a particle moving in a quartic confining 
potential. 

We remark that work on the ABC model on the ring has been carried 
out recently by Bodineau et al. [8]. They investigated analytically the pair 
correlations, and thus also the local density fluctuations, in the uniform 
state, and showed that the latter diverge as T \ Tc. They also suggested 
a possible formula for the large deviation and discussed the possibility of 
deriving this from the macroscopic fluctuation theory [9J. The model on 
the ring was also studied, in a different context, by Fayolle and Purtlehner 
[10) . who obtained results which agree, for the case in which the system on 
the ring has the same stationary Gibbs measure as does the system on the 
interval, with those obtained here. Results concerning the relation between 
the hydrodynamic equations for the ABC model and the LDF were obtained 
by Bertini et al. [9]. 

The outline of the rest of the paper is as follows. In Section [2] we de- 
fine the Gibbs measure with respect to which the dynamics satisfy detailed 
balance. We then study some properties of this measure; in particular, 
we describe its ground states, which are nonunique whenever two species 
have equal mean densities which are greater than or equal to 1/3. We in- 
vestigate also some properties of the microscopic correlation functions at 
finite-temperature: the mean field nature of the interactions in this system 
leads in the iV — > oo limit to local measures which are exchangeable. 

In Section [3] we consider the scaling limit of the model and briefly dis- 
cuss the nature of possible limiting density profiles. In Section |4] we ob- 
tain the Helmholtz free energy .7^ as a functional of the scaled densities 
Paix), X G [0,1] (this is equivalent to obtaining the LDF) as well as the 
Euler-Lagrange equations satisfied by minimizers of the proof that mini- 
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mizers exist and must satisfy these equations is postponed until Section [TOl 
In Section [5] we investigate the minimizers of T by studying in detail the 
solutions of the Euler-Lagrange equations. These solutions, which describe 
all stationary points of T ^ are given by periodic (elliptic) functions describ- 
ing the motion of a particle in a one dimensional quartic potential. There 
are many such solutions for large /?; despite this, we prove uniqueness of the 
globally minimizing density profiles for all /3 so long as the mean densities 
are not all equal. Details of the uniqueness proof are given in Sections [6] and 
El In Section [8] we discuss the phase diagram of the model and describe a 
perturbation expansion around the uniform state. In Section [9] we show that 
is a convex functional of the /9a(x)'s (with fixed r^'s) at high temperature, 
which gives an alternate proof of the uniqueness of minimizers of T in this 
regime. A concluding discussion is given in Section [TTl 



2 The steady state of the model 

A configuration C, of the ABC model on the interval is an A^-tuple 
(Ci, • • • , Civ), with = A., B, or C. We will let r]a{i) be a random variable 
which specifies whether a particle of species a is present at site i: r]a{i) = 1 
if Ci = a and r]a{i) = otherwise, so that r]A{i) + ??b(«) + ??c(*) = 1 and 
Na = ^j^iVaii)- The time invariant measure for the weakly asymmetric 
dynamics described in Section [1] is a canonical Gibbs measure, 

iypiO = Z~'exp[-(3EM{0], (2-1) 

where 

^ N-l N 

^^^0 = nY.Y. ivc{i)VB{j) + VA{i)vc{j) + VB{i)VA{j)] (2.2) 

i=l j=i+l 

and Z is the usual equilibrium normalization factor, that is, the (canonical) 
partition function with fixed particle numbers N^- To verify this invariance 
one checks detailed balance for the dynamics: if C is a configuration with 
particles of types a and a + 1 on sites i and i + 1, respectively (recall the 
convention that a, a + 1, ... run cyclically through A, B, and C), and 
is the configuration with these particles interchanged, then the transition 
rates are that ( a,t rate e~^/^ and ^ C at rate 1, and the 

detailed balance condition e~^f^ -vpi^C,) = l-vp{C,'^'^^'^) follows from ()2.ip and 
(|2.2p . For [3 < oo, i^p is the unique stationary measure; this follows from the 
transitivity of the dynamics on the set of all Nl/ (NaINbWcI) configurations 
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consistent with YliVa{i) = ^a- When /? = all these configurations are 
equally likely. 

One can rewrite (j2.2p in various forms which diff'er from each other only 
by functions of Na, Nb, and Nc, and this does not affect vp when the Na 
are fixed. As an example, an equivalent energy is 

. N N~l N ^. 

^*^ = n\ Ej'[^^(^')-^a(j)]+3 5] VB{i)VA{j)\. (2.3) 

j=l i=l j=i+l ^ 

We shall generally use the form E]\f{Q given in (j2.2p . as it clearly exhibits 
the cyclic symmetry between the different species. 

Remark 2.1. (a) The energy of (j2.3p is of particular interest when 
the B species is not present, i.e., when Nb = 0. In this case we may 
call the A particles just particles and the C particles holes, and the model 
reduces to the weakly asymmetric simple exclusion process (WASEP). Then 
E"^ = —N^^ jr]A{j), which is the energy arising from an external field of 
magnitude pushing the particles to the right. See [llj and Section 7.2.2 
of the review [12] for discussions of the partially asymmetric simple exclusion 
process on an interval. 

(b) On the ring, the stationary measure for the WASEP gives equal weight to 
all (jv^) configurations. In contrast to the situation on the interval, however, 
the dynamics is not reversible, i.e., does not satisfy detailed balance, with 
respect to this stationary measure. We thus have a true nonequilibrium 
stationary state (NESS). 

2.1 Ground states 

For (3 = oo there are many configurations stationary for the dynamics: 
all those of the form • • • ABCAB • • • , where a boldface letter denotes a 
block of particles of the corresponding species. These configurations are 
local minimizers of the energy i?jv(C), in the sense that any interchange of 
adjacent particles of different species will raise the energy, but not all of them 
are obtained as limits of when /? ^ oo. In this limit the measure will be 
concentrated on the ground states of -Eat, that is, on the global minimizers 
of E]y{() [13]. As we will see, the energy per particle in the ground state C 
is eo = N~^En{C) = min{r^rB,r^rc,rBrc}, where = N^/N. 

We describe the ground states by considering several cases. It is helpful 
to note that (12. 2p may be summarized as saying that there is a contribution 
of to the energy each time a B particle lies to the left of an A particle. 



7 



a C to the left of a i?, or an ^ to the left of a C We also remark that if 
a local minimizer contains a sequence BCAB of four blocks with a total of 
ka particles of type a, then one may lower the energy by regrouping these 
particles into three blocks, BCA, ABC, or CAB, unless fc^ = fee > ks 
(and similarly for ABC A and CABC). 

(i) If one species is not present, say = 0, then the ground state energy 
is zero and there is a unique ground state configuration with all particles 
of type a + 1 to the left of all particles of type a + 2. This is the WASEP 
discussed in Remark 12.11 In this case there are no local minimizers other 
than the global minimizer. 

From now on we suppose that none of the vanish. 

(ii) If one of the A^^^ is greater than the other two, say > maxjiV^, Nc} 
(other cases can then be found by cyclic permutation), then it is easy to see 
that there is a unique ground state C consisting of three blocks ordered as 
ABC, with eo = rArc- 

(iii) If two of the Na are equal and the third is smaller, say Na = Nc > Nb, 
the ground state will be Nb + 1 fold degenerate, consisting of three or four 
blocks ordered as BCA, CAB, or BCAB, that is, some of S's will be at 
the left side of the interval, followed by all the C"s, then all the ^'s, and 
finally the remaining i?'s. Here eo = rATB = tbtc- 

(iv) If all of the Na are equal, Na = Nb = -/Vc, then the ground states are 
the shifts, by an arbitrary number of sites, of the configuration with three 
blocks ordered as ABC; see the discussion in Section [1] of the correspondence 
in this case with the model on the ring. The ground state is thus A'^-fold 
degenerate, and eo = 1/9. 

The ground state phase diagram of the system in the right triangle where 
ta^tb ^ and va + tb < 1 is given in Figure [J (remember that rc = 
l — rA — VB)- The ground state is unique everywhere except on the three line 
segments originating from r^i = = 1/3 and terminating at the midpoints 
of the sides of the triangle, corresponding to cases (iii) and (iv) above. On 
the edges of the triangle the system reduces to the WASEP. 

2.2 Finite temperature correlations and local states 

The correlation functions at finite (3 are 

)), (2.4) 

where (•) denotes the expectation with respect to the Gibbs measure vp. We 
point out here some simple relations that these satisfy, both in finite volume 
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Figure 1: The {rA,rB) ground state phase diagram. It exhibits three phases 
separated by three first order hues. In region (a), the A density is larger 
than that of the B's and C's and hence the order is such that the A-block is 
in the middle. Similarly (b) and (c) correspond to regions where B and C 
densities are largest. 



and in the limit N ^ oo, Na/N r^. (With some abuse of notation we will 
use ra to denote both the density N^/N in a finite system and the limiting 
value of the density in the N ^ oo limit.) 

First, it is easy to check from (|2.ip and (|2.2p that, whenever all the Na 
are positive, 

and this implies that for all A^, 

3 3 

ll{r]aiN)) = ll{va{l)). (2.6) 

a=l a=l 

Note that when r = (1/3, 1/3, 1/3) and the system on the interval coincides 
with that on the ring, ()2.6p is a consequence of the translation invariance of 
the latter. 

Next we note that there is a constant K, which depends only on /? and 
the Va, such that for fixed x G (0, 1), ji and j2 integers with < ji < j2, 
and ao, ai, and a2 distinct particle species, 

lim {r,^,{[xN\)l^^,{[xN\ +n)r]^,{[xN\ + h)) = K, (2.7) 

N—^co 
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where [xA^J denotes the largest integer less or equal to xN. To see that 
the limit in (j2.7p is independent of x, ji, j2, ao, qi, and 02 we first write 
io = [xN\, ii = \xN\ + ji, and 12 = \xN\ +32-, and observe that if C, 
is a configuration with C,i^ = ao, Cn = oti-, ai^d = 02, and ^ is the 
configuration on — 3 sites which is obtained from C, by omitting the sites 
ii, 12, and ^3 and renumbering, then En{C) = En-3{S,) + 1 + A; here the 
terms 1 + A arise from those terms in the energy (I2.2p involving the particles 
on sites ii, 12, and ^3, and A, which depends on ^, the i^, and the a^, satisfies 
|A| < j2/N. Thus from I^J^ and ([22]), 

Z-'Mio)VaAh)VaA^2)) = ^e-^^--©+i+^. (2.8) 

The corresponding expression for a different choice of the ik and Ok will 
differ from (j2.8p only in the replacement of A by some A', from which the 
result follows. 

The results (|2.6p and (j2.7p are closely related: each shows that a product 
of three occupation numbers, or their expectations, is site independent. We 
will see another such relation in Section[5) for all scaled density profiles Pa{x) 
which are stationary points of the free energy functional (see Section |4]) the 
product pa{x) pb{x) pc{x) is in fact independent of x. 

Similar identities can help us understand the nature of local states [H] 
in this model. By the local state p(^x) position x G (0, 1) we mean the 
infinite volume state, if it exists, whose correlation functions for particles of 
type ai, . . . ,afi at sites [xN\ + ji , . . . , [xN\ + are given by 

h 

\im^(jlri^^{lxN\+j,)). (2.9) 

i=l 

(If such a local state does not exist then one may identify subsequences 
Nk along which all the limits (j2.9p do exist, and thus obtain a family of 
local states.) An argument similar to that establishing (12. 7p shows that for 
X £ (0, 1), ji, . . . ,jh fixed, and P any permutation of {1, ... , h}, 

(UllVaAlxN\+j,)) 

lim V = 1- (2-10) 

This implies [15] that //(^j) is a product measure or a superposition of product 
measures, which of course reflects the mean field nature of the interaction. 
Dynamically this corresponds to the fact that when — > 00, the exchanges 
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become symmetric and the invariant measures are just superpositions of 
extremal product measures [11]. (In fact this last argument apphes also on 
the ring with unequal densities r^.) 

3 The scaling limit 

We now turn to the properties of the model in the macroscopic scaling limit 

N ^oo, No^/N r„, i/N x. (3.1) 

We say that there is a unique (macroscopic) density profile in this limit when 
there exists a function p{x) = {pa{x), psix), pci^)) such that [16j 



lim Prob 

N-*oo 



{x)ipaix) dx 



><5^=0, (3.2) 



for all a, all > 0, and all piecewise smooth functions ipa{x) on the interval 
[0,1]. Such Pa{x) must, of course, satisfy the conditions 

< /9a(x) < 1 and / pa{x)dx = Va for a = A,B,C, 

•^0 (3.3) 
^/)c.(x) = 1. 

a 

Equation (|3.2p expresses the convergence, in a certain sense, of the ran- 
dom profile A^~^ r]a{i)di/]\[ to the non-random profile Pa{x)'- specifically, 
for each the random variable '?a(^)Va(^/-^) converges in the sense 

of (|3.2p to the non-random quantity / pa{x)ipa{x) dx. It may happen, on 
the other hand, that the limiting profile is itself random, in the sense that 
there is some family of limiting profiles p{u}) = {/9q(ca;;x)}, indexed by the 
variable to lying in some sample space Q and each satisfying (13. Sh . and a 
measure k on 17, such that for each test function ipa{x), 



Pa{-]x)(pa{x)dx (3.4) 



in the sense of convergence in distribution. (In fact, (j3.2p is also convergence 
in distribution, to a non-random limit.) We would then say that the limiting 
profile is not unique and that the system can exist in more than one phase, 
i.e., is in a phase transition region. This is exactly what happens to the 
ferromagnetic Ising model in two or more dimensions in a periodic box at 
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temperatures below the critical temperatures and fixed magnetization m £ 
(— m*,m*), where m* is the spontaneous magnetization [T7j. In our system 
this happens, as we shall see, if and only if = = rc and T < (27r\/3)~^. 
More generally, it may be that the limit in (j3.4p will not exist; in that 
case, one would expect to obtain a limit in the given form by passing to a 
subsequence. 

4 The free energy 

The question of whether or not there exists a limiting profile(s) as given 
by (13. 2p or (13. 4p . and of which of these alternatives applies, is related to 
the question of the existence and uniqueness of minimizers of the LDF for 
this system [IS]. If n{x) = {na{x))a=A,B,c denotes a general continuum 
density profile satisfying (|3.3p (written in terms of n rather than p) then the 
logarithm of the probability of finding the profile n{x) in the scaling limit 
is asymptotically —NTi,B{{n{x)}, where .Fld is the LDF. Possible limiting 
density profiles are thus those which minimize ^ld- 

To obtain from the microscopic measure the probability of finding the 
profile 77.(2;) it is necessary (roughly speaking) to sum 1^/3 (C), as given in (j2.ip . 
over all (microscopic) configurations C consistent with n(x); this summation 
will yield a restricted partition function divided by the full partition function 
Z. Taking the logarithm of this ratio, dividing by N, and taking the scaling 
limit of the result will then yield 

:FLBiMx)}) = HMx)}) - mi T{{n{x)}), (4.1) 

n(x) 

where P~^J-{{n{x)}) is the Helmholtz free energy of the system restricted 
to having a density profile n{x): 

HMx)}) = mMx)}) - S{{n{x)}). (4.2) 

Here the macroscopic energy £{{n{x)}) and entropy S{{n{x)}) are the limits 
of the microscopic energy and entropy per site; as indicated, these limits 
depend only on the profile n{x). For the energy, this is due to the mean 
field nature of the microscopic energy En of ()2.2p . For the entropy, which 
is the logarithm of the number of microscopic configurations consistent with 
n{x), this is due to the Bernoulli nature of each component of the local 
measure, which implies that all microscopic configurations consistent with 
n{x) have the same energy and hence the same weight. Thus 

•^({"D = /5 / dx dz[nA{x)nc{z) +nB{x)nA{z) +nc{x)nB{z)] 

Jo Jx 
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The profiles n{x) minimizing (j4.3|) — that is, the candidates for hmiting den- 
sity profiles — will represent a compromise between the entropy, which wants 
to keep all the densities uniform, and the energy, which wants to keep the 
different particle types segregated. 

The next theorem establishes the existence and properties of these min- 
imizers. 

Theorem 4.1. Suppose that < /? < oo. Then given positive numbers r^, 
tb, and rc with rA + fB + fc = 1, let F{rA,rB,rc) denote the infimum the 
values of J^{{n}) over (measurable) profiles n{x) satisfying the constraints 
(j3.3p . Then there exist infinitely differentiable functions paIx), Pb{x) and 
pc{x), possibly non-unique, that satisfy these constraints and achieve the 
minimum: 



Moreover, there is a 6 > so that each pa satisfies 5 < pa{x) < I — 5 for all 
x; thus the minimizer is an interior point with respect to the constraint of 
taking values in [0, 1] and so satisfies the Euler- Lagrange equations (ELE) 
obtained from T . 

The proof of this theorem is given in Section [TOj We note here, how- 
ever, that the minimizers lie in the interior of the constraint region essen- 
tially because the entropy term has infinite normal derivative at the bound- 
ary. For /3 = oo, when only the energy counts, the densities minimizing 
lim^^oo (3^^ J' = £ will be the continuum limit of the ground state configu- 
rations described in section 12.11 and need not satisfy the ELE. 

To obtain the Euler-Lagrange equations satisfied by the minimizing pro- 
files we must take the variational derivatives of J- {{n}) with respect to two 
of the density profiles while maintaining = 1 and na{x) dx = ra', 

here we will treat n^(x) and nB{x) as independent, with nc{x) = l—na{x) — 
nB{x). Defining !Fa{x) = 5!F / 5na{x) to be the variational derivative taken 
as if the profiles n^(x), nB{x), and nc{x) were independent, we will then 
have 



^{{p}) = F{rA,rB,rc) ■ 



(4.4) 



61^ 



Tb - Tc- (4.5) 



bnA 



nc=i-nA-nB 



Sub 



nc=l-nA-nB 
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Imposing the constraints na{x) dx = leads to the ELE 

— = constant, J^b — J^c = constant, (4.6) 

where 

Ta{x) = logna{x) + P [na+i{z) -na+2{z)]dz + 1 + Pra+2- (4.7) 

Jo 

Simple manipulations then show that under the constraint '^^^^^i-^) ~ ^ 
the derivatives J^a satisfy 

^»^^o. = 0, (4.8) 

a=A,B,C 

which implies, rather surprisingly, that solutions of (|4.6p will satisfy J-a{x) = 
constant for all a, i.e., that the functional derivatives of J- can be taken as 
if the ria were independent. From (14. 7p . then, the minimizing profile p{x) 
must satisfy the ELE 

p^{x) = p„(0)e'5/o'[p-+2(?/)-P"+i(?/)l'^y, (4.9) 

for all X £ [0, 1]. It follows that 

Pa(l)=/>a(0)e^('-"+^-^'^+^), (4.10) 

which is consistent with (12. 5p . The ELE may also be written in differential 
form: 

'^P'' -PpAipc-pB), (4.11a) 



dx 
dps 

dx 
dpc 

dx 



f3pB {pA - Pc) , (4.11b) 
ppc {pB - pa) ■ (4.11c) 



These equations were derived in [7] for the case r = (1/3,1/3,1/3). One 
obtains minimizing profiles by solving these equations with the constraints 

1 

Pa{x) dx = Tq,. (4.12) 







14 



Remark 4.2. As indicated above, one may consider as a function of 
UA and ub alone, with nc{x) = 1 — ^^(a;) — nsix). The corresponding 
variational derivatives Aa(2;) and Ab(x) are the local chemical potentials for 
the density profiles, 



Xb = — ^c- 



(4.13) 



which are constant at a stationary point of J- (see (j4.6p ). These are given 
by 



5nA{x) 



Ayi({"'yi}, {ns}) 
nA{x) 



log 



6nB{x) 



nc{x) 
XsiinA]: {nB}) 
nB{x) 



13 {l-2,nB{z))dz + l3{rc-rB), (4.14) 



log 



(3 {l-2,nA{z))dz + (3{rA-rc), (4.15) 



nc{x) 

with the replacements nc{x) = 1 — nA{x) — nB{x) and rc = 1 — — r^. 



5 Solutions of the ELE equations 

We now turn to a detailed study of the solutions of the Euler-Lagrange equa- 
tions (j4.1ip with the constraints ()4.12p . As already noted in Theorem 14. H 
for any given {3 and r^, t^, there must exist at least one solution of these, 
corresponding to a minimizer oi T. There may, however, be many solutions; 
these will all correspond to stationary points of and more than one of 
these stationary points may be a minimizer. In this section we determine 
all the minimizers. 

5.1 The WASEP 

We first consider briefly the case in which one of the densities is zero, say 
= 0, when the system reduces to the WASEP (see Remark 12. ip . This 
system, which can be solved explicitly on the interval, has in the scaling limit 
a unique limiting density profile — i.e., minimizer of the free energy — given 
by the solution of (14. lip with pb{x) = 0: 



Pa{x) = ^^^^ = 1 - pc{x), xG[0, 1]. 



(5.1) 
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The constant D is determined by the constraint pA{x)dx = r^i = 1 — rc- 
When /? — > cxD, (j5.ip reduces to the ground state configuration in which all 
the A particles are pushed to the right: 

, Jo, ifO<x< 1-r^, 

°° 11, if 1 - ryi < X < 1. 

The local measure //(^.^ in this model is, as expected, the Bernoulli measure 
with density pa{x)- (On the ring the stationary state of the WASEP is a 
nonequilibrium one in which all configurations have equal weight and the 
limiting scaled density profile is constant.) 

5.2 Properties of solutions. 

Through the remainder of this section we suppose that all the are strictly 
positive. Let us begin by discussing the properties of some given solution 
p{x) = [pa{x)., pb{x), pc{x)) of the ELE equations (j4.1ip with the con- 
straints (14.121) . By differentiating \og{pA{x) pb{x) pc{x)) with respect to x 
we see that the ELE imply that there is a constant K such that 

Pa{x) pb{x) pc{x) = K, < X < 1, (5.3) 

and that the condition 

Pa{x) + pb{x) + pc{x) = I, 0<x<l, (5.4) 



is preserved by the equations. Equation (15. 3p is the scaling limit version of 
(12771). 

Since pa{x) and pc{x) have sum 1 — pb{x) and product K/ pb{x) they 
must be the two roots of the equation + (1 — pB{x))r + K/pb{x) = 0, 
from which 



f f \ ^ pb{x){1 - pb{xW - "^K 

Pa{x) - pc{x) = ±\ -—^ . (5.5) 

V Pb{x) 

Then squaring ()4.11bp and using (j5.5p we find that Pb{x) is a solution of 

p'{xf + ^P^Uk{p{x))=Q, (5.6) 

where 

UK{p)^\Kp-\p\l-pf. (5.7) 
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Of course, the same argument shows that pAi^) and pc{x) must also be 
solutions of (|5.6p . 

To study (j5.6p we write t = 2j3x and let y{t) = p(t/2(3); then y satisfies 

^y'itf + UKiyit))=0. (5.8) 

This is the equation of the zero energy solution of a mass 1 particle moving 
in a potential Uk- Because Uk is quartic in y the solutions are elliptic 
functions; see Appendix |Al For K > 1/27, Uxiv) is strictly positive for 
< y < 1 and so no solutions with y in this range exist. For K = 1/27, 
Uxiy) has a local minimum with value at y = 1/3, and (j5.8p has constant 
solution y{t) = 1/3. For < K < 1/27, Ur has four zeros, 0, a{K), h{K), 
and c{K), where 0<a<6<l<c and UK{y) < for a < y < 6. 
See Figure EJ Since we are interested in solutions of ()4.1ip which satisfy 
< Pa{x) < 1 we consider only the solutions of (j5.8p which oscillate between 
a and h. Let yx denote the solution of (|5.8p which satisfies yxi^) = «; yxit) 
then has period 

= 2 r —^== . (5.9) 

L ^/-2UK{y) 

For < t < yxit) is determined by inverting the relation 

, , fy dw , , 

t = iiy) = / ; (5-10) 

Ja ^/-2Uk{w) 

yxit) is then obtained for all t by extending to an even function of period 

TK- 

Remark 5.1. We note some properties of this solution; those listed in (a) 
are elementary, and we relegate the proofs of (b) and (c) to Appendix IbI 

(a) yxit) is even, is periodic with minimal period tk, has local minima at 
integer multiples of tk and local maxima at half integer multiples of tk , and 
is monotonic between these points. Moreover, yxit) = yx^s) if and only if 
either s + t or s — t is an integer multiple of tk- 

(b) lim;^ /'i/27'''K = 47r\/3 = 2/3c and limi^\o(''"x/ ln(l/-fC)) = 6. Moreover, 
for any e > 0, 

^ r dy , , 

lim — T-TT / , = 2, 5.11 

Mln(l/i^)A y'-2UK{y) 



and 



lim — ^ r = 1. (5.12) 

K^o ln(l/K) ^-2UK{y) 
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-0.01-' 



Figure 2: Plots of Uxiv) for K = 1/20, 1/27, 1/50, 1/100, and 1/1000. 



This means that in the hmit K ^ oo the trajectory i/Kit) will spend two- 
thirds of its time near the origin and one-third near y = 1; see for example 
the third trajectory shown in Figure HI 

(c) tk is a strictly monotonic decreasing function of K for < < 1 /27. 
Now for a = A,B,C there must be a phase shift ta such that 

Pa{x)=yK{2p{x-l/2)+ta), < X < 1, (5.13) 

that is, each pa{x) is obtained by looking at the solution yK{t) within a 
window of length 2/3 centered at some value ta, and rescaling from t to x. 
The phase shifts are not independent; in fact, 

tA = tB + TK/3 and tc = tB-TK/3 (5-14) 

(see (27) of |7]). In verifying (I5.14p the requirement (I4.12p of certain average 
densities is not relevant. What matters is that we consider three solutions of 
the ELE (j4.1ip or equivalently (with the rescaling t = 2f3x) three solutions 
yK,a{t) of satisfying 

v'ka = VkAvk^c - yK,B)/2, 

v'k.b = VkMvka - yK,c)/2, (5.15) 

yK,c = yK,c{yK,B - yK,A)/2, 

with yK,a{t) = yxit+ta)- We may assume that these are defined for all t and 
without loss of generality that ts = 0. It is helpful to view the trajectories 



18 



in the y-y' phase plane; see Figure [3l At time t = the B trajectory hes at 
the point on the y axis marked Bl. Since the velocity on the B trajectory 
is zero at this point, it follows from (j5.15p that yK,A{0) = yK,c{0), that 
is, the A and C trajectories he at points such as those marked Al and CI, 
respectively; again from (jS.lSh we know that Al must be in the upper half 
of the phase plane, since y'x a(^) ^ow follow the motion for a time 

period At until yK,A reaches the y axis at point A2; since now y^^(At) = 
it follows that the B and C trajectories are at points B2 and C2 with equal 
y coordinates, and it follows immediately (using the symmetry around the y 
axis) that the travel time along each of the six orbital segments in the figure 
is At. Thus At = t^ = t^ + 2 At = ts + tk /S, and tc = ts — tk /'i- 



y' 





^^B2 




— ■ 

Al 






\a2 




\0.1 0.2 


0.3 0.4 


0.5 


0.6 


0.7 


J 0.8 






V 
















.CI 









Figure 3: Phase plane orbit of yxit) for K = 1/100. The time intervals 
between the six marked points are all equal. 



Let us introduce the notation 

{z)i = ^J^z{t)dt, (5.16) 

where \I\ = Jj dt, for the average of the function z over the set / (which 
will always be a union of intervals), and define Yxit) = {yK)[t-i3,t+/3]- Then 
with (l5T3l) and (f5l4l) . KT2t\ becomes 

YK{tB) = rB, YK{tB+TKm=rA, YK{tB-TK/5)=rc. (5.17) 

The problem of solving ()4.1ip and ()4.12p is now the problem of finding K 
and tB satisfying (j5.17p . If either K = 1/27 so that yxii) = 1/3 is constant, 
or 2(3 is an integer multiple of tk, then Yxit) = 1/3 is constant and (I5.17P 
has a solution (with arbitrary ts) if and only if = ?"_b = rc = 1/3. 
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Figure m shows the curves yxit) and yxit ^tk/S) for several values of 
K. To obtain a solution of (j5.17p for some K one views the corresponding 
three curves in a window of length 2/3 (corresponding to the full lattice, i.e., 
to the original unit interval after the variable change in (j5.13p ) centered at 
tB- In this context we label a solution by an integer which is one more than 
the number of full periods (plus, perhaps, a fraction of a period) fit into the 
window: we say that a solution p{x) of ()5.17p with K < 1/27 is of type n, 
for n = 1,2,..., if 

in-l)TK <2P <nTK. (5.18) 

We do not assign a type to the constant solution yx = 1/3, which exists for 
K = 1/27, VA = tb = rc = 1/3, and every value of /3, as discussed above. 

5.3 Uniqueness of solutions 

The next theorem summarizes our results about solutions of the ELE and 
answers the questions posed in the first paragraph of this section. We con- 
sider the problem of finding yx and t b satisfying (j5.17p . 

Theorem 5.2. (a) If ta = rB = rc = 1/3 then there exist (i) the constant 
solution, (a) for (3 > n[3c = 27rn\/3, n = 1, 2, . . a solution of type n unique 
up to translation, and (Hi) no other solutions. The minimizer of the free 
energy is, for (3 < [3c, the (unique) constant solution and, for (3 > f3c, the 
type 1 solution. 

(b) For values of r other than (1/3, 1/3, 1/3) there exists for all (3 a unique 
type 1 solution which is a minimizer of the free energy. 

In the remainder of this section we prove part (a) of this theorem and 
give an overview of some parts of the proof of part (b). More technical 
parts of the proof of (b) are given in Sections [6] and [71 Before beginning, we 
summarize some simple properties of the function Yk of (15.171) . 

Remark 5.3. (a) As remarked above, if either K = 1/27 or 2(3 is a multiple 
of TK then Yxit) = 1/3 for all t. 

(b) For all K,(3 other than those of (a) a plot Y^it) and Y^it it ri^/3) will 
look much like one of the sets of curves in Figure HI In particular, since yK 
is even and r/^'-periodic, so is Yk, and from Y^{t) = {2l3)~^[yK{t + /3) — 
yK{t — (3)] it follows that Yk is strictly monotonic between integer multiples 
of tk/2. From this we see that if YK{t) = Yk{s) then either s + t or s — t 
must be an integer multiple of tk- 

(c) To find all triples (r^, rB,rc) for which (j5.17p has a solution for a given 
K it suffices to consider values oi tB satisfying < < r/^/G. For such tB 
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O.Sn 




Figure 4: Plots of yxit) (solid), yK(t + (dotted), and yxit — tk/^) 

(dashed) for K = 1/50 (top), K = 1/6400 (middle), and K = 5-^2'^° = 
1/27487790694400 (bottom). 



the possible triples satisfy rs < rc < ?'a; other orderings of the same sets 
of values are found for other ranges oi ts- 

Proof of Theorem \5.S\( a): When r = (1/3, 1/3, 1/3), (|5.17p becomes 

YxitB) = YK{tB + tk/3) = YkHb - tk/3) = 1/3. (5.19) 

For every /3, one solution of ()5.19p is the constant VK{t){= Pa{x)) = 1/3 
corresponding to = 1/27. Moreover, if 27rn\/3 < f3 < 27r(n + l)\/3 then 
by Remark 15.11 there exist n one-parameter families of nonconstant solutions 
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(indexed hy ts) obtained by solving 

2l3 = jTK,, j = l...,n, (5.20) 

for Kj. Note that the solution (j5.20p is of type j. No other solutions are 
possible, because (I5.19P is inconsistent with the final observation in Re- 
mark [5i3][^b). 

To complete the proof we note that, as is easily seen, perturbing the 
uniform solution via Pa(x) — > 1/3 + e cos(27r(2; + Xq,)), where xa = xb + 1/3 
and xc = Xb — 1/3, decreases the free energy when /3 > /3c, so that this 
solution cannot be a minimizer, and that we will prove in Section [7] (see 
Theorem I7.ip that no solution of type n with n > 2 can minimize the free 
energy. □ 

Proof of Theorem \5.S^f h) : It was shown in Theorem 14 . 1 1 that for any positive 
r^, tbi and rc which satisfy Yla ^« = 1 there exists at least one minimizer of 
the free energy and that this minimizer satisfies the ELE. By Theorem l7.ll 
no solution of type n with n > 2 can minimize the free energy, so this 
minimizer must be of type 1. It remains to prove uniqueness of this solution; 
we will give the full proof in Section [6] (see Theorem 16. ip . but here we 
illustrate the idea of the proof by sketching the argument for the special case 
tb < 1/3, = rc > 1/3, in which symmetry considerations considerably 
simplify the discussion. 

We suppose then that for some such r and some /3 there are two distinct 
type 1 solutions of (I5.17p . and derive a contradiction. It follows from rA = rc 
and Remark I5.3l fb) that the phase shift tB must be zero for each solution, 
so that we are looking at some zi = yxi and Z2 = yK2 on the interval 
J = [-j3,f3] which satisfy (see (f5l^ ) 

{zi)j = {Z2)j = rB. (5.21) 

We take K2 < Ki and write 9i = ri^./2, so that (because these are type 1 
solutions) (3 < 61 < 02- The situation is as in Figure 15.31 the qualitative 
features of this figure, which we use below, are obtained from ^2(0) < 2^1(0) 
(which holds because y_ft:(0) = a{K), the smallest positive zero of the poten- 
tial Uk{p) (|5.7p . decreases with K), the symmetry, and the fact that (|5.2ip 
implies that the curves zi and Z2 must cross in the interval J. In the full 
proof below we will verify that the curves cross exactly twice, as shown. 
We write Ii = [— 6*1,6*1] and I2 = [—^2,^2] and observe that 



{zi)h = {Z2)h = 1/3- 



(5.22) 
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Since {z2)i2 is a (weighted) average of {z2)j = rs < 1/3 and {z2)j2\j, 
necessarily 

{Z2)l,\j > TB- (5.23) 

Moreover, 

(^2>/2\j > {Z2)h\j > {zi)i,\j. (5.24) 

where the first inequality holds because Z2 is decreasing to the left of the 
origin and increasing to the right, and the second because Z2 > zi on Ii \ J . 
But then 

= -^rB + —^{zi)i,\j 

= rB+ {{zi)i^\j -rB) 



< rB + —^ {{z2)h\j-rB) 



< r-B + {{Z2)i2\j - ^b) 

= {Z2)i„ (5.25) 
which contradicts (j5.22p . □ 
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6 Uniqueness of type 1 solutions 

Our goal in this section is to prove 

Theorem 6.1. For any /? > and any positive ta, rs, and rc which 
satisfy Yla = ^ and are not all equal to 1 /3 there exists at most one 
type 1 solution of (jS.Sp satisfying (|5.14p - (j5.17p . 



Throughout the section we will consider two solutions yKi , yK2 of ([57 
with < K2 < Ki < 1/27; we write Tj rather than tk^, i = 1,2, for the 
corresponding periods, and note that ti < T2 by Remark 15. li fe). We begin 
with a preliminary result. 

Lemma 6.2. (a) For any to G M and n G Z the curves yxiit + ^o) and 
yK'A't) intersect exactly once in the interval (nr2/2, (n + l)r2/2). 
(h) (i) IfO<t< T2/6 then 

yK,{2T2/3 -t)> yxA'^n/S - t), (6.1) 

and (a) if in addition T2 < 3ri then 

yK2{T2-t) <yKi{Ti-t). (6.2) 

Proof, (a) Suppose that n is even so that yK2 is increasing on the interval 
(nr2/2, (n + l)r2/2); the proof for n odd is similar. For all t and to the zeros 
a{K) and b{K) of Ur (see ([SZD) satisfy 

a{K2) < a{Ki) < yKA^ + t^) < b{Ki) < b{K2), (6.3) 

and since yK2{f^'T2/2) = a(i^2) and yK2{{^+'^)'^2/2) = h{K2), the existence of 
at least one intersection follows from the intermediate value theorem. On the 
other hand, if yKiit + to) = yK2{t) = 2/ foi' some t G (nr2/2, (n+l)r2/2) then 
from dsn) and \y'j^^{t + to)\ = ^-2UKAy) < V^2U^^ = y^,(i), 

and this is inconsistent with the existence of two such intersection points. 

(b) It follows from (a) that for any ti,t2, yKi{t + ti) and yK2(t + ^2) can 
intersect at most once on any interval of monotonicity of yK2{t + *2)- Thus 
it suffices to verify that (j6.ip and (|6.2p hold for t = and t = r2/6. For 
(|6T]) this follows from yK2{ri/2) = h{K2) (see ([63])) and 

/2r2X l-a(i^2) , l-a(i^i) /2rn 

yK2 [-^) = — ^ — > — ^ — = yK. [-^) . (6.4) 

For (j6.2p it follows from yx2(''"2) = a{K2) and 

/5t2\ 1 - KK2) ^ 1 - b{Ki) /^5ri\ / T2\ 

yK2 [-^) = ^ < ^ = yK. [-^) < yK. [n - ^) , (6.5) 
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where at the last step we have used monotonicity of yxi on [ri/2,ri] and 
the condition T2 < 3ti. □ 

In the remainder of this section we suppose that, given appropriate phase 
shifts, the solutions yxi and considered above provide two type 1 so- 
lutions satisfying (j5.17p - ()5.14p for the given r^, rs, and rc, and from this 
derive a contradiction. It is convenient to view these solutions in the same 
interval, which we take to be J = [— /?, This means that there are phase 
shifts ti and t2 such that if Zi{t) = yx^it + U) and Wi(t) = yx^it + + Tj/3) 
for i = 1,2, then 

{zi)j = {z2)j = TB and {wi)j = {w2)j = va- (6.6) 

By Remark 15. 3f c) we may again assume that > > rs and that 

0<U<^, i = l,2. (6.7) 
b 

The fact that these are type 1 solutions, and our assumption that Ki > K2, 
imply that 

The key idea in the proof is the same as that for the simple case consid- 
ered in Section [5.31 but the general case presents two additional difficulties. 
First, without guidance from symmetry the choice of the intervals Ii and I2 
is more delicate, and in fact we must consider two distinct cases in which 
they are chosen by different prescriptions. Second, one must in some cases 
apply the reasoning to zi and Z2, and in others to wi and 11)2, and one must 
show that one or the other is possible. Sorting out these cases requires some 
detailed analysis of the geometry of the curves. 

Proof of Theorem \ 6.1[ Equation ()6.6p implies that the curves zi (t) and 
Z2{t) must intersect at least once in the interior of the interval J, and since 
the length 2/3 of J is less than T2, Lemma 16.2( a) implies that they cannot 
intersect more than three times. On the other hand, from (j4.10p it follows 
that zi{(5) / Z2{j3) = zi{— 13) / Z2{— 13) . It is convenient then to consider two 
possible alternatives: 

zl: Z2{I3) > zi{(3), Z2{—(3) > zi{—(3), and the curves zi and Z2 
intersect precisely twice within J; 

z2: Z2{(3) < zi{(3), Z2{—(3) < zi{—(3), and either both inequalities 
are strict and the curves intersect precisely twice within J, or both 
are equalities and the curves have one intersection within J. 
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A similar analysis leads to corresponding alternatives for wi and W2- 

wl: W2{I3) < wi{P), W2{—(3) < 'Wi{—(3), and the curves wi and ^2 
intersect precisely twice within J; 

w2: W2{P) > W2{—(3) > tfi (—/?), and either both inequalities 

are strict and the curves intersect precisely twice within J, or both 
are equalities and the curves have one intersection within J. 

We will first show that if zl (respectively wl) occurs then the argument of 
Section [5^ may be applied to zi and Z2 (respectively wi and ^2) to derive 
a contradiction. Then we show that either zl or wl or must occur. 

Now observe that necessarily f3 > t2- For otherwise, < t2 — P < t2+ f3 < 
T2/3, so that Z2{t) = yK2{t + ^2) is increasing on J and by Lemma [6.2r a) 
cannot intersect zi{t) more than once there. It is then helpful to further 
subdivide the situation into two cases in which we can explicitly identify 
intervals on which Z2 is monotonic: 

Case 1: f3 + t2> T2/2. In this case we write J = Ji U J2 U J3, where 

Ji = [-13, -t2], J2 = [-t2,T2/2-t2], and J3 = [T2/2 - t2, 13]. 

Case 2: (3 + t2 < T2/2. Here we write J = Ji U J2, where Ji = 

[-(3,-t2] and J2 = [-t2,(3]. 
Then Z2 is decreasing on Ji and J3 and increasing on J2. 

Suppose now that zl occurs. Since Z2{—I3) > zi{—P) and Z2(— ^2) = 
yK2{^) < zi{—t2) it follows from Lemma lOT a) that zi{t) and Z2{t) intersect 
once on Ji. We now consider separately the two cases introduced above. 

Case 1: Define h = [P - n,/?] and h = [P - T2,/?], so that J C h C h- 
Since /j has length Tj, we again have ()5.22p and therefore (15.230 . Moreover, 
()5.24p follows from the monotonicity of Z2 on l2\ J and then the fact that 
Z2 > z\ on 1\ \ J, which follows from this monotonicity. Lemma 16.2( a). and 
the fact that Z2 intersects z\ on J\ . But then we may deduce a contradiction 
between (16. 6p and (15.221) . just as we did in (15.251) . 

Case 2: The argument is similar. Let I2 = [—t2 — T2/2, T2/2 — 12] = I2, 
where /| = [-^2 - T2/2, -t2] and = [-^2,^2/2 - t2]. Let h = ![ U If, 
where l[ = [—a, — ^2] and /J = ["^2, b], be an interval of length ri satisfying 
J C /i C /2, with a and b chosen so that 



Now by arguing as in Case 1 we find that (15.241) holds with I2 and Ii replaced 
either by I2 and or by I2 and /[. Averaging these two equations (the 
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weights for this averaging, 



are independent of i by (|6.9p ) we again find that (j5.24p itself holds, and the 
argument proceeds to a contradiction as in Case 1. 

If wl occurs then we may derive a contradiction similarly: we need only 
replace Zj by Wi throughout and change the sign of some of the inequalities. 

It remains to show that either zl or wl must occur; to do so we assume 
that z2 occurs and show that then wl must also. The occurrence of z2 
implies that we must be in Case 1, i.e., that 

/3 + t2>j, (6.11) 

since if Z2{—P) < zi{—P) there can be no intersection of zi and Z2 on Ji and 
hence the two intersections of these curves must occur on J2 and J3, while 
if 2:2 = Zi{—P) then there are three intersections which must occur on 
Ji, J2 and J3. We claim further that necessarily 

- ^ > - (6.12) 
For if (|6.12p does not hold then 

z2iP)=yK,{P + t2)>yK, (/^ + *i + ^^^^^-^) >yK,{P + ti) = zi{p). 

(6.13) 

Here for the first inequality we have used ()6.7p , ()6.8p , (|6.1ip , and the falsity 
of (I6.12p . which together imply that 

f </3 + t2</? + ti + ^^^<f , (6.14) 

and then the monotonicity of on [T2/2, 2r2/3]. For the second inequality 
we use Lemma l6.2r b.i). which is applicable by (j6.14p . But now 

W2iP) = yK2{f3 + t2 + J 

< VK, ( /? + i2 + n - ^ 

< yXi(/3 + ti+ri/3) =u;i(/3), (6.15) 
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which imphes wl. Here the first inequahty is from Lemma l6.2( b.ii). apph- 
cable because ()6.7p . ()6.8p . and (|6.1ip imply 

< 2r2/3 - /3 - t2 < r2/6, (6.16) 

and (j6.1ip implies that /3 > T^j'i and so T2 < 3ri/2, and the second from 
(j6.12p and the monotonicity of yxi on [2ri/3, ri], which is applicable by 
(|6.16p and the observation that (|6.1ip implies that (3 > T2/3 > ti/3. □ 



7 Solutions of type n, n > 2 

It is easy to see that type n solutions of (j5.17p with n > 2 do exist for 
certain but not all values of /3 and r. The next result, however, shows that 
these solutions are not of physical interest. 

Theorem 7.1. No type n solution of (j5.17p with n > 2 minimizes the free 
energy. 

We begin by rewriting the free energy in terms of the variable t = 2(3x. 
Consider then triples z{t) = {zA{t), zs^t), zc{t)) of functions, defined in 
some interval [c, d] of length 2/3 and satisfying < Za{t) < 1 and ^„ Za{t) = 
1. The free energy T{{z{t)}) of z is (see ()4.3p ) 

Him) = m-H£{{z{t)}) - s{{zitm, (7.1) 

where 

dt ds [zA{t)zc{s) + ZB{t)zA{s) + zc{t)zB{s)] (7.2) 

is the energy and 

S{{z}) = -2 dt [zA\ogZA + zslogZB + zclogzc] (7.3) 

the entropy. We will use a rearrangement procedure for these triples. 

Let c = t{) < ti < • • • < tra = d he a. partition of [c, d\ and for k = 
1, . . . , m let /fc = (tfc_i, tfc]- Given a permutation a of {1, 2, . . . , m} we may 
rearrange the intervals Ik in [c, d\ , together with the restrictions of z to 
each interval, into the order Ia[i) i • • • i Ia{m) ■, thus defining in the obvious 
way a new function w on [c, d], the rearrangement of z. A formal definition 
of w (which we will not use in the sequel) may be given as follows: The 
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permuted intervals arise from the partition c = sq < si < ■ ■ ■ < Sm = d, 
with Sj = Yl{i\a{i)<j}i^i — ti-i) , and this permutation is implemented by the 
map V' : [c, d] — > [c, d] defined by 

ipic) = c, ipix) = X - tk-i + Sk-i if X £ Ik. (7.4) 

The rearrangement of z is then w = z o tp^'^, i.e., Wa{'^{t)) = Za{t). 
The entropy of z and w are the same. The energy of z is 

m 

1=1 i<*<i<'Ti 

where 

£i= I dt [ ds[zAit)zcis) + ZBit)zAis) + zcit)zBis)], (7.6) 
Jteii J seii,s>t 

and for i < j, with ZjQ, = -Za(i) dt, 

Sij = dt ds[zAit)zc{s) + ZBit)zA{s) + zc{t)zB{s)] (7.7) 

= Zi^AZj,C + Zi^BZj^A + Zi,CZj,B- 

It is natural to think of (j7.7p as expressing iS({2:(t)}) as a sum of "self 
energies" of the restrictions of z to the intervals and "interaction energies" 
between the portions of z in different intervals. The energy £{{w{t)}) is 
obtained similarly; the self energy contribution will be the same and the 
interaction energies will differ only for interval pairs whose order is 

interchanged by the rearrangement; thus 

:F{{w{t)}) - :F{{z{t)}) = £{{w{t)})-£{{zm 

E A£,,j, (7.8) 

{i,j\i<j, o-(j)<(T(i)} 

where 

A£ij = [Zi,CZj,A + Zi,AZj,B + Zi,BZj,C - Zi,AZj,C - Zi,BZj,A - Zi^cZj,B]- (7.9) 

Proof of Theorem \7.1\ Suppose that the triple z{t) as above gives a solution 
of (j5.17p on the interval [c, d]; specifically, this means that for some K and 
some phase shift ts, c = tB — (3, d = tB + P, and ZA{t) = yxit + tk /S), 
ZB{t) = yxii), and zc{t) = y_ft:(t — t/^/3). Suppose also that the solution 
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is of type n for some n > 2, which imphes in particular that t = tk < 2/?. 
We will show that for some rearrangement w of z the change in free energy 
()7.8p is negative. Since clearly {wa)[c,d] = {za)[c,d] ^or all a, w will satisfy 
the same constraints as z; thus z cannot be a minimizer of the free energy 
under those constraints. 

The desired rearrangement will be defined in terms of a partition 
to, . . . ,tm with to = c, tk = (ko + k)T/3 for A; = 1, . . . , m — 1, and tm = d; 
here ko is the largest integer such that fcoT/3 < c and m is the smallest 
integer such that {ko + m)T/3 > d. Note that here Ii and Im have length at 
most t/3 and Ij, j = 2, . . . ,m — 1, has length exactly t/3. It is easy to see 
that (m - 2)t/3 < 2/? < mr/3, so that from (fSTTSj) . 3n-2<m<3n + l 
and in particular m > 4. 

Now for each j, j = 1, . . . , m, there will be an a G {A, B, C} such that 
Za{t) > 2Q±i(t) for all t £ Ij (see Figure we will then say that Ij is a full 
a-interval if it has length t/3 and a partial a-interval if it has length less 
than t/3 (which is possible only if j = 1 or m). The types of the intervals 
are in cyclic order (again see Figure S]): if Ij is an a-interval then /j+i is an 
(q + l)-interval. Recalling that m > 4 we focus on four consecutive intervals 
Ii, . . . , Ii+3 . Without loss of generality we may assume that /j is an A- 
interval, and we will then write Ii = Ai, /j+i = B, /j_|_2 = C, and /j+a = Ar, 
i.e., the initial configuration is AiBCAr- Consider then the permutation of 
these intervals in which we first exchange B and C, then interchange Ai 
with C (they are now adjacent) and A^ with B, and finally exchange Ai and 
Ar, leading to C,Ar,Ai,B. We will show that this rearrangement lowers the 
energy, completing the proof of the theorem. 

The calculation is simplest when all four of the intervals are full, so we 
consider that special case first. For a full a-interval Ij we have Zj^a = C 
Zj,a±i = C, where 

rT/3 /"'"/3 /"'"/3 

(= ZA{t)dt, and C = / ZB{t)dt= / zc{t)dt, (7.10) 

and so C > C- We may then compute from (jT.Op and (jT.lOp the energy 
difference AiSjj arising from reversing the order of a full Oj-interval Ii and 
a full ttj -interval Ij, where i < j: if aj = ai then A£ij = 0, while if 
aj = ai ± 1 then ASij = ±e, where e = {( — C)^ > 0. Rearranging the 
intervals as above thus first increases the energy by e, then decreases it by 
2e; the final exchange of Ai and Ar leaves the energy unchanged. The net 
change of energy is — e. 

We now consider the general case in which either Ai or Ar or both may 
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be partial intervals. After a translation, if necessary, of [c, d] by a multiple 
of r we may assume that Ai = (u, t/3], B = {t/3,2t/3], C = (2r/3,r], 
and Ar = (r, r + v\ for some u,v with {) < u < v < t the condition 
u < V comes from the fact that if A^ is full then u = and if Ai is full then 
V = r/3, while if both are partial then v + t — u = d — c = 2P>t. The 
energy change from the rearrangement is the sum, with j = z + 3, of 

• — C)^ from the exchange of C and B, 

• — (C — C)i^i,A — ^i,B) from the exchange of A^ and C, 

• — (C — C)i^j,A — ^j,c) from the exchange of B and Ar, and 

• Zi^cZj,A + Zi,AZj,B + Zi^BZj,c - Zi,A^j,C - Zi,BZj,A - Zi,cZj,B from the 
exchange of Ai and Ar. 

This energy may conveniently be expressed in terms of 

V^a^ = Za{t)dt, 7/^2)=/ z^{t) dt, and r]^^^= Za{t)dt, 



(7.11) 

z,^^ = ^^^^\ c = r?\^^+^r + ^r. (7.12) 



usmg 

,(2) , „(3) ^ . _ „(1) , J2) . _ iX) , „(2) (3) 



and 

C = r/g^ + r?(,'^ + = 7?g) + r?g) + (7.13) 

(see Figure [6]) . Finally, if we introduce the strictly positive quantities 9^^^ = 
^ — ryg^ and = r]'^ —rl~^\ then the total energy change may be written, 
using the identity displayed in (j7.13p . in the manifestly negative form 

□ 



Remark 7.2. Energy considerations for full intervals correspond closely 
to those of single particles: cyclic order ABC is energetically favored, the 
energy cost of a reversed pair BA, CB, or AC is independent of the species 
involved, etc. The conclusion above that CAAB is favored over ABC A 
parallels the conclusion of Section 12.11 that the ground state of a particle 
system with a majority of A^s has a block containing all the A's between 
blocks of the Cs and the 5's. For partial intervals, energy considerations 
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Figure 6: Plots of za, zb, and zc in the interval [0,Ti<-/3] for K = 0.005, 
with typical values of u and v. 



are more subtle; the final step in the rearrangement above, the interchange 
of Ai and A^, is irrelevant when these intervals are full but is needed when 
they are partial because Ar is richer in species C than in B and the reverse 
is true for Ai. 



8 Phase diagram of the ABC model 

Our results about the minimizers of for specified P, r^, and rs provide 
the following picture of the canonical phase diagram of the model (see Fig. 
[7]). The three first order lines in the T = plane, discussed in Section El do 
not extend to finite temperatures. There is a first order line for = = 
rc = 1/3 which starts at T = and terminates at the critical temperature 
i.e., as one crosses this line by varying the overall densities 
r^, the density profiles change discontinuously for T < Tc. On the line of 
equal densities the transition at Tc is second order; the divergence of the 
variance in local density fluctuations as one approaches Tc from above on 
this line was investigated by Bodineau et al. via a study of the 1/A^ 
dependence of the correlation function. 

To see how the density profiles change as a function of the average den- 
sities near the = 1/3 line we take = rc and vary rs- The density 
profiles change continuously at high temperatures /? < /3c as the profiles 
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Figure 7: The (T, r^jr^) phase diagram: The temperature is added as a 
vertical axis to the right triangle. The T = phase transition lines do not 
extend to T 7^ 0, except at the equal density point = 1/3 where it extends 
as a first-order line at low T and ends at the critical point f3 = 2tt\/3. 

become uniform at Tq, = 1/3. However, when /3 > (3c, the B particles prefer 
to be in the middle section of the interval for rs > 1/3 and to be symmet- 
rically arranged at both ends for rs < 1/3, and this change of profile is 
discontinuous. To see this discontinuity we define an order parameter, 

/■3/4 

VTB = 2 / pB{x)dx, (8.1) 

Jl/4 

and calculate it in the steady state. Figure [8l shows the order parameter 
ttb as a function of for /3 = 10, 12 and 15. Clearly, for low temperatures, 
(3 > (3c, TTs shows a discontinuity at density vb = 1/3. 

8.1 Perturbation expansion 

In this section we apply a perturbation expansion to calculate the density 
profiles for average densities which deviate slightly from 1/3. In particular 
we carry out a perturbation expansion of the profiles for average densities 
rA = 1/3 + € + 7] , vb = 1/3 — 2e and rc = 1/3 + e — 7/ with small parameters 
e and rj . For e = rj = the density profiles are homogeneous (namely, 
Pa{x) = 1/3) at high temperatures, and they become non-homogeneous at 
(3 < (3c = 27r\/3 . We now expand the density profiles around the homoge- 
neous densities: 

PAix) = 1 + a{x), pb{x) = 1 + b{x), and PC (x) = I + c{x) (8.2) 
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Figure 8: The order parameter ttb defined in (j8.ip is shown as function 
of e = {tb — l/3)/2 for different temperatures /3 = 10,12 and 15. Other 
parameters are, = rc = (1 — rB)/2. Data were obtained via a mean-field 
dynamical solution in a system of 200 sites. 



with small deviations a(x), b{x) and c{x) which satisfy a{x) + b{x) + c{x) = 0. 
As a result of this constraint, one is left with two independent functions 
which we choose to be b{x) and s{x) = a{x) — c{x). The average densities 
are then fixed by, 



/•1/2 ri/2 

/ b{x)dx = — 2e and / s{x)dx = 2r], (8-3) 

J-l/2 J-1/2 



where, for convenience, the interval is taken to be —1/2 < x < 1/2. In terms 
of b{x) and s(x) the equations for the minimizers (j4.1ip reduce to 

-T- = -s + (3bs 
ax 6 

We proceed by expanding b{x) and s{x) in terms of two small parameters 
u and f , which, at the end of the calculation, will be determined by e and 
r/ . To be explicit, let us write 

oo oo 

Kx) = ^K{x), and s{x) = ^Sfc(x), (8.5) 

fc=l k=l 
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where bk{x) and Sfc(x) are functions of x which are of order k in the smah 
parameters u and v. To first order in u and v ()8.4p yields the fohowing 
equations for bi{x) and si{x) 

dbi 13 dsi 

— = -si and — = -f3bi . 

ax 6 ax 

Solving these equations one finds 

bi{x) = ncos(Qx) + usin(ax) 

si{x) = —uV3sm{ax) + vV3cos{ax) (8-6) 

with a = P/y/S . Expanding b{x) and s{x) to second order in u and v the 
following equations are obtained for 62(2;) and ^2(3;), 

db2 1 3a 2 2\ / \ 

= -j=as2 —{u —V ) sin(2ax) + 3an?; cos(2qx) 

— ^ = —V^ab2 H — ^ (u^ — v^) cos(2qx) + 3\/3q;to sin(2ax) .(8.7) 
dx 2 

These equations can be solved to yield 

b2{x) = 2^^^ ~ cos(2a2;) + uvs\n{2ax) 

S2{x) = —^{u^ — v'^)si\i{2ax) — \/2>uv cos{2ax). (8.8) 

The equations for the third order terms, 63(x) and S3(x) , are then 

db^. 1 3 2 ■ / ^ 3 2 ( \ 

dx v3 2 2 

= — \/3a63 + '^^ anw^ cos(aj;) + avuP' sin(ax) , (8.9) 
ax 2 2 

where w'^ = + v^. To proceed with the analysis one eliminates S3(x) 
from the equation for 63 (x) to obtain the following second order differential 
equation for 63 (x) 

—r-T + oP'^'i = 'ioPuw'^ cos(ax) + 'ia^vvP' sin(ax) . (8.10) 
ax^ 

This equation may be solved to yield 

.,3 9 ,,3 o 
03(x) = -auw xsin(ax) — -avw xcos(qx) . (8.11) 
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Inserting this solution in (|8.9p results in the following solution for S3 



S3{x) = ^ aunP'x cos(aa:) H — ^^^avnP' s\ii{ax) . 
In summary, to third order in u the density profiles are given by 



^.12) 



h{x) 



s{x) 



ucos{ax) + f sin(ax) + ^(^^ ~ '^^) cos(2ax) + sin(2ax) 

3 2-/\^ 2 / \ 

+ -auw xsui[ax) — -avw xcos[ax), 
—uVSsm^ax) + vV3cos{ax) H — ^(^^ ~ sin(2aa;) 

— \/3nt> cos(2ax) H — ^auw^x cos(ax) 

3\/3 2 / \ 
H —avw xsin(ax). (8.13) 



To third order in u these equations may be reexpressed as 
b{x) = ucos(7j;) + u sin(7j;) H — (n^ — v^) cos(27x) 



s[x) 



+ tiv sin(27x) + 0(ti ,v ) , 



V3, 



— n\/3sin(7x) + vV3cos{jx) + ~^(^^ ~ sin(27a;) 
VSuv cos(27x) + 0(n^ w"^) , (8.14) 



with 



j = a--aw 



.15) 



Note that s(x) = 6(x + 27r/37) — 6(x — 27r/37), which is consistent with 

The expansion parameters u and t; may be expressed in terms of e and 
r] using (|8.3|) . This leads to 



2e 



27/ 



-sin(7/2) 

7 





" 1 


n + 


— sm7 
.27 . 



(tx2 - v^) 



7 



sin(7/2) 





'V3 . 


V — 


sm 7 

7 



lit; . 



^.16) 



For e = r/ = a nonzero solution of (j8.16p exists only if sin 7/2 = 0, and 
thus by (j8.15p only if a > 2n, i.e., /? > Pc- For /? < f3c, then, all profiles are 
homogeneous. For (3 > j3c there exist type n solutions with 7 = 2n7r and 
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w small whenever (3 is just above n/?c- Their amplitude w = (u^ + w^)^/^ 
is determined by /? via (|8.15p but u and v are otherwise undetermined, 
corresponding to the translation invariance of the set of solutions. 

For non-vanishing e or r/ it is convenient to expand (j8.16p to third order: 



2e 



2r/ 



+ 



+ 



— sin(a/2) 
a 



u + 



1 

2a 



■ sm a 



■ cos(a/2) H — sin(a/2) 
a 



uw 



2\/3 

a 

3^/3 



sin(a/2) 



^/3 



■ sm a 



a 



-cos(q/2) + -sin(Q/2) 
a 



uv 



vw 



.17) 



Then (j8.17|) has one or more non- vanishing solutions for u or w at any tem- 
perature. The equilibrium (minimizing) solution is the one which is of type 
1. The minimizing density profiles, as has been shown, are non-homogeneous 
both above and below f3c] they vary continuously with the temperature /3, 
and no phase transition takes place. 

The small u and v expansion may be used to verify the first order nature 
of the transition at e = = just below the critical point. To demonstrate 
this point we take, for simplicity, = and consider small e and a = 27r+Aa 
with Aa > 0. In this case u = 0, the amplitude u satisfies 



27r 



and the density profile is given by 



2e = — (Aa - 37ru^)n 



(8.18) 



h{x) = u cos(27r -|- Aa — 37rn^ 
At e = 0, equation (j8.18p has three solutions, 



(8.19) 



u 







Aa 
37 



(8.20) 



Since the stable solution satisfies Aa — "ii^v? < 0, it follows from (j8.18p that 
for the stable solution one has e/u < 0. Thus for e > the negative u 

- Y^Aa/37r, while for e < the positive u 



e\0 



U 



solution is stable and lim, 
solution is stable and lim, 
discontinuous at e = 0; see Figure [H 



Aa/37r. The density profile is therefore 
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9 Convexity of the free energy at high tempera- 
ture 

It follows from Theorem 15.21 that if /3 < /3c = (27r\/3) then for any positive 
rA, TB and rc with = 1 there is a unique solution of the ELE (I4.11|) 

with the constraints (j4.12p : that is, a unique stationary point of the free 
energy T under those constraints. One would then like to know also the 
nature of the fluctuations about this equilibrium profile, and in particular 
to show that these fiuctuations are Gaussian with a finite covariance. To 
establish this it is necessary to show that the free energy has a positive 
definite second variation at the minimizer in the space of density profiles 
satisfying the above constraints, which implies that it is a strictly convex 
functional in a neighborhood of the minimizer. We believe that this is 
indeed true generally, with the obvious exception of the line segment P > Pc, 
T~A = '''B=T~c = ^/'^- In this section we show a global convexity property 
in a restricted temperature range : for ^ < 47r/3 = (4/3\/3)/3c, T (with the 
mean density of species a constrained to be r^) is globally strictly convex. 

To establish the convexity we compute the second variation. Consider 
the free energy functional J-{{n{x)}) on the convex set of nonnegative den- 
sities satisfying 

nA{x) + nsix) + nc{x) = 1 and / na{x)(ix = ra- (9-1) 

Jo 

To make the variations, let (pA and (/i^ be bounded continuous functions 
with ^ ^ 

(l)A{x)dLX = / (l)B{x)dx = . (9.2) 



/O JO 

Then for small values of t, 



(nA, nB,nc) (ua + t(j)A, ns + t(j)B,nc - t[(t)A + 4>b\) (9.3) 

is an admissible variation, and all admissible variations are of this form. 
The order contribution from the entropy is 



- i dx 

2 Jo \_nA[x) 



-^4>a{x) + —^(t>l{x) + —^{<t^A{x) + ct)B{x)f 
nA[x) nB{x) nc[x) 



> 



Jo 



' dx[cl>\{x)+^lix)] , (9.4) 
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where we have used the fact that for any positive a, b, and c with a+b+c = 1 
the matrix 



(9.5) 



is positive semidefinite, since its diagonal entries are nonnegative and its 
determinant is {abc)~^[{a — 6)^(1 — c) + c(l — c/2)^ + 3c'^/4]. Because the 
energy is quadratic, the order contribution from the energy is independent 
of {nA,nB,nc); it is 



(3 dx dye(y-x)[-0A(x)(0A(y) + 0B(y)) 
Jo Jo 

+ ^B{x)(l)A{y) - (M^) + 0B(x))</.B(y)] . (9.6) 

However, since for functions / and g with Jq f{x)dx = Jq g{x)dx = 0, 

"1 pi pi pi 

dx / dyQ{y - x)f{x)g{y) = - i dx I dy Q{y - x)g{x)f{y) , (9.7) 
) Jo Jo Jo 



reduces to 



3/3 / dx / dyQ{y - x)(l)Bix)(t)Aiy) 



{91 



Now let TC be the Hilbert space of square integrable functions / on [0, 1] 
with f{x)dx = 0, and define on operator K onTi 



Kf{x) = P (^1^ 



f{y)<iy 



(9.9) 



where P is the orthogonal projection onto Ti in L^([0, 1]). Then it is easy 
to see that 



dx / dy@{y - x)(l)Bix)<j)Aiy) = {4>b, K(t)A)'H ■ 



(9.10) 



Combining the entropy bound (j9.4p with (j9.8p and (|9.10p . we see that the 
second variation is bounded below by 



dx{<t)\{x) + <t)l{x)) + W{<t)B,K<t)A)H ■ 



(9.11) 



We need to show that this quadratic form is nonnegative definite for (3 < 
47r/3, which we shall do with a spectral calculation. 
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To do the calculation, define 
V^nix) = \/2 sin(n27rx) and 
Then one has 



ipn{x) = \/2 cos(n27rx) . 



K(pn{x) 



1 



n27r 



and 



n2TT 



If we now write 



^(anV^n + bnlpn) and (/)B = ^(C„(^„ + dnlpr, 



n=l 



n=l 



then ()9.1ip becomes 



E 

n=l 



3/? 



(a„ + 6„ + c„ + d„) + — -(-c„6„ + d„a„) 



n27r 



(9.12) 
(9.13) 

(9.14) 
(9.15) 



This is nonnegative as long as /3 < 47r/3, which proves the strict positivity 
of the second variation, and hence the convexity of J-' for such f5 

We also remark that if in the equal density case ta = tb = tq = \ 
one considers quadratic variations around the constant profile then one can 
improve the lower bound (j9.4p on the second variation of the entropy to 



(9.16) 



3 / {4>A + 4>B + (t^A<i>B) dx. 



Some modification of the remainder of the argument then shows that the 
constant profile is indeed a local minimum all the way to f3c- 



10 Existence of minimizers 

In this section we prove Theorem 14. H that is, we show that there exist 
profiles p{x) which minimize the free energy functional J-'{{n{x)}) of ()4.3p 
and that these satisfy the corresponding Euler Lagrange equations. The 
latter cannot be taken for granted, as there may be situations in which 
the minimizing density profile {p(x)} is on the boundary of the permissible 
domain, e.g., one might have PAix) = for some values of x, and when 
that happens the minimizing profile need not satisfy the ELE. In the proof 
of existence we will consider the weak topology on profiles, in which a 
sequence {fk} converges to /, where /fc, / € L^([0, 1]), if and only if 

lim / dx fk{x)(j){x) = I dxf{x)(t>{x) (10-1) 

k-~*'X'Jo Jo 
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for all bounded and measurable functions (j). We will use the fact that, by 
the Dunford-Pettis Theorem [19], the subset /C := {/ | < f{x) < 1 a.e.} of 
-L^([0, 1]) is compact in this topology (since it is clearly uniformly integrable). 
We need also one other preliminary result. 

Lemma 10.1. The free energy T is lower semicontinuous in the weak 
topology, that is, if n = {nA,nB,nc) and Uk = {nk,A,nk,B,nk,c), 
k = 1,2,..., are profiles with each sequence {n^^} converging to Ua in 
the sense of (110. ip then liminffc^oo •^({'i-fc}) > -^({^D- 



Proof. It suffices to show that if lim„_ 
weak topology then 



,fn = f and lim^^oo 5n = in the 



lim 

n— >oo 



dx / dzQ{z - x)fn{x)gn{z) 



dx 



dzQ{z - x)f{x)g{z) 
(10.2) 



and 



lim inf 



dx f n{x) In fn{x) > / dx f{x)lnf{x) 



(10.3) 



To verify (jl0.2p . define hn{x) = Jq dzQ{z — x)fn{x) and h{x) = dzQ{z 
x)f{x), and notice that by the definition of weak convergence. 



lim hn{x) = h{x) 



:io.4) 



almost everywhere. Then as each /i„ takes its values in [0, 1], the dominated 
convergence theorem implies that /i„ converges to h strongly in the norm. 
Fix e > 0, and pick N so that \\hn — h\\i < e for all n> N. Then 



< 



dxgn{x)hn{x 
1 

dx {gn{x) - g{x))h{x) 



/ dx g{x)h{x 
Jo 

+ 



dx gn{x){hn{x) - h{x)) 



. (10.5) 



But 



dx gn{x){hn{x) - h{x)) 



< Ibnilooll/in - h\\i < e since < gn < 



1. Likewise, since < h < 1, liuin^oo Jq dx {gn{x) — g{x))h{x) = 
3y the definition of weak convergence. Thus for all n sufficiently large, 
dx gn{x)hn{x) — dx g{x)h{x) < 2e, and (jl0.2p follows. 
Equation (jl0.3p is an easy consequence of the convexity of the function 
—S{f) = dx f (x) In f (x) on K. Indeed, for any real number a, the set 
K-a = {f ^ : —S{f) < a} is convex. It is also closed in the norm 
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topology, for if {/n} is a norm convergent sequence in fCa with limit / 
then there is a subsequence converging almost everywhere, and passing to 
the subsequence and using the dominated convergence theorem, with the 
obvious bound \ fn\ < 1/e, yields —S{f) = liuin^oo — S (fn) < a, so that 
/ G ICa as well. But by a theorem of Mazur [20], a convex set in a Banach 
space is closed if and only if it is weakly closed, and so ICa is weakly closed. 
It follows easily from this that if {/„} converges weakly to /, then 

- S{f) < liminf -S{fn) , (10.6) 

n— >oo 

which is (fTOS]) . □ 



Proof of Theorem \4-.l\ Let {n^.^, fc, nc^^} be a minimizing sequence for 
.F, with rifc^Q G ^! ~ ^"^^ So ^o(^) = ^cn so that 

lim J^{{nk{x)}) = F{rA,rB,rc) (10.7) 

fc— >oo 

By the weak compactness of /C we can choose a subsequence along which 
each no,,k converges to some pa G fC. From (jlO.ip it follows immediately that 
Jq dx Pa{x) = ra, from na,k = 1 that Yla Pa = ^, and from Lemma flO. II 
that 

^{{pix)}) < liminf.F({nfc(x)}) = F{rA,rB,rc) ■ (10.8) 

Thus, J^{{p{x)}) = F{rA,rB,rc). 

We next show that each pa is uniformly bounded below by some constant 
(5 > 0. Suppose for example that this is not the case for pA, so that for each 
(5 > 0, the set = {x\ pa{x) < 6} has strictly positive measure \Ds\. We 
can then "fill in the hole" in pA- we transfer some mass from the larger oi pB 
and PC to in a way that preserves the constraint pAix)+ pb{x)+ pc{x) = 
1, and then transfer mass in the opposite direction on a certain "safe" set M 
to restore the constraints p^dx = 1. Because "entropy abhors a vacuum" , 
these transfers will, for 5 sufficiently small, strictly lower the value of thus 
contradicting the optimality of {/>«}. 

On the set M, pA should be strictly bounded below; we take M = {x \ 
Pa{x) > and note that |M| > because pA has mean value r^. Next, 

let Es be any subset of Ds with < \Es\ < |M|, let Fg be the subset of Es 
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on which ps > pc-, and let Gs = Es\Fg. Finally, define {pa} by 

Pa{^) = Pa{x)+S1e,{x)-6^-^^1m{x) , (10.9) 

Pb{x) = pb{x)-51f,{x) + 8^-0^1m{x) , (10.10) 

pc{x) = /9c(x)-51g,(x)+5^1m(x) . (10.11) 

Using the fact that, on D^, max{/9^(x) pb{x)} > (1 — 5)/2, and on M, 
Pb{x), pc{x) < 1 — pa/2 < 1, one sees easily that if < 5 < min{rA/2, 1/3} 
then Pa satisfies < ^^(x) < 1. 

It is now easy to verify that for all 5 sufficiently small, 

- S{pa) = -S{pa) + {SH6) + 0{6))\Es\. (10.12) 

Moreover, changes in S{pb), S{pc), and the energy components of are 
of order 0{5)\Eg\, and thus, for 5 sufficiently small, J^{{pa}) < J^{{pa}), 
which contradicts the optimality of {pa}- Hence for some 5 > 0, it must be 
the case that Dg is a null set. The same argument applies of course to pB 
and Pc, and hence each of the {pa} is uniformly bounded below by a strictly 
positive constant. Then from the fact that pa{x)+ pb{x) + pc{x) = 1 almost 
everywhere, it follows that Pa{x) < 1 — 25 for almost every x and each a. 

Finally, since we have an interior minimum, {pa} must satisfy the ELE 
for J^, that is, J^a — and Tb — must be constant (see ()4.6p ). which 
from (j4.7p implies that pa/ Pc and Pb/pc are differentiable. From this the 
constraint ^aPaix) = 1 implies readily that all the pa are differentiable, 
and the argument of Section [3] them shows that they satisfy the ELE in the 
form (|4.9p . This implies infinite differentiability. □ 



11 Discussion and Conclusion 

The dynamics defining the weakly asymmetric ABC model on the interval 
are entirely local and identical to those for the model on a ring, except 
for the insertion of a barrier between sites and 1. Remarkably, when 
Na = Nb = Nc = N/3 this has no effect on the stationary state of the 
system, which (as was known before for the ring) is given by a Gibbs measure 
with mean field type asymmetric long range pair interactions (see (|2.ip . 
(j2.2p ). When the species numbers are unequal, however, the stationary 
state on the ring is a nonequilibrium one, with a net current, while the state 
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on the interval is again given by a canonical Gibbs measure with the same 
interactions. 

In the scaling limit — > oo, N^/N — > r^, the system on the ring 
with TA = tb = Tc = ^/2> was known to have a second order transition 
ai (3 = Pc = 27r\/3 from a macroscopic state with uniform density profiles 
Pa{x) = 1/3 to one in which the density profiles are periodic, with period 
one. This transition, which by the isomorphism between the ring and in- 
terval systems at equal densities carries over to the interval, corresponds to 
a broken symmetry in which the phase of the typical density profile in the 
scaling limit is uniformly distributed over [0,1]. This corresponds in (|3.4p 
to = [0, 1], with the measure k uniform on Q. 

In the case of the ring there is strong evidence [8] , but no proof so far, that 
there will be a transition from uniform to nonuniform typical configurations 
also for a range of unequal, but strictly positive, densities. It might then be 
natural to conjecture that something similar happens on the interval. Our 
results here, however, prove that this is not the case: for unequal densities 
the change from /3 = to /? = oo is smooth, with unique density profiles 
Pa{x) for each /3 < oo. In the limit /3 — > cxd the system will segregate into 
three or four blocks of different species; see the discussion in Section 2.1 of 
the ground states of the finite system. When < < ta = '''c the limiting 
configuration has four blocks of particles, BCAB, with the B particles 
evenly divided between the left and right ends of the interval, that is, the 
degeneracy of the ground state of the finite system, in which the B particles 
could be arbitrarily divided between the two ends, is broken in the scaling 
limit when (3 ^ oo. 

A striking characteristic of the resulting (T, r^,rB) phase diagram 
(Fig. 7) is that one observes a phase transition at precisely one value of 
the external parameters (here ta and tb) as the temperature is lowered. 
Such a feature is typically encountered for grand canonical ensembles of 
systems in the presence of an ordering field (e.g. the (T, Hx, Hy) phase dia- 
gram of the XY model in a magnetic field H = {Hx, Hy)). Note, however, 
that here we are considering a mean field model and permitting only those 
changes in the system in which the densities r^i, r^, and rc are kept strictly 
fixed. The question of what happens for more general variations, that is, if 
we consider a grand canonical ensemble, has to be investigated separately 
(see note 5 below) . This is unlike the case of systems with short range inter- 
actions, where one has an equivalence of ensembles. What is interesting in 
this model compared with the usual mean field situation is the rich spatial 
structure of the equilibrium states, which is due entirely to the directional 
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asymmetry of the interactions. Thus if we modified the energy function 
in (|2.2p to involve a sum over all j ^ i, not just j > i, we would obtain 
= NaNb + NaNc + NbNc which, in the canonical ensemble, would 
give uniform minimizing density profiles at all temperatures, since the in- 
teraction would not depend on the spatial structure and the entropy term 
prefers the uniform state. 

As we have shown, the minimizing profiles are given by pieces of the 
graphs of elliptic functions uk describing the r-periodic trajectory of a par- 
ticle with zero energy in a quartic confining potential. It follows from the 
general analysis of the Euler-Lagrange equations that these elliptic func- 
tions have the property that yK(t) + yxit + t/3) + yKit — t/3) = 1 while 
yK{t)yK{t + T/3)yK{t — t/3) = K is independent of t. It is not clear to us 
whether or not this property of these elliptic functions yxii) is known in 
the literature. 

We note that, as pointed out in Section O the local measures must 
be either product measures or superpositions of product measures. Clearly 
when ()3.2p holds, as it will wherever the minimizers are unique, = ^p{x) 
will be a product measure with density while for r^i = = rc = 1/3 
and f3 > Pc the Hx will be a superposition of the measures iJ,p{z) as z varies 
uniformly over [0, 1]. 

We end this section with several open problems. 

1. It follows from our analysis that for large (3 there will be solutions of the 
ELE of type n > 1 which are stationary points, but not global minimizers, 
of We have not determined whether or not these correspond to local 
minima. 

2. It follows from the proof (see Section[9]) of the strict convexity of J^{n{x)} 
for small /3 that the fluctuations about the minimizing p(x) at fixed are 
(constrained) Gaussian. We know that this is not true for (3 > f3c and 
T~A = T~B = rc = 1/3, due to the existence of a one parameter family of 
minimizing profiles. What of the fluctuations about minimizing densities 
for other values of when (3 is not very small? 

3. As mentioned in the introduction, Evans et al. considered the ABC model 
on the ring with general rates g^^^ (q ^ 7) and found that the stationary 
state is Gibbsian whenever A^~-^ log(q'o+i,a+2/^a+2,a+i) is independent of a. 
This is possible, of course, only when the ratio of any two of these logarithms 
is rational, and then only when is a multiple of some smallest possible 
system size, just as for the model considered in the body of this paper the 
stationary state can be Gibbsian only if is divisible by 3. To study the 
weakly asymmetric scaling limit one would take qa,a+i = <la+i,a 
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with some normalization, say '^a = 1- The energy 

i=l j=i+l a 

would then be translation invariant when Na = CaN. As in the current pa- 
per, however, one may use (jll.ip to obtain a Gibbs measure exp[— /JE'at] 
on the interval for any values of the N^, and then study the scaling limit 
()3.ip . It would be interesting to carry out an analysis of the free energy 
minimizing profiles also for this model. 

4. The ABC model can be extended to higher dimensions either on a torus 
or in a box [21j . There does not, however, seem to be any natural way to 
extend the Gibbs measure to these systems for any values of the r^'s. 

5. Our analysis of J-'{{n}) has been carried out entirely for the canonical 
ensemble, the natural one for the exchange dynamics. From the point of 
view of the Gibbs measure ()2.ip it is natural to consider also the grand 
canonical ensemble in which the are not fixed. In that case one would 
have to consider variations in the grand canonical partition function, the 
pressure, for fixed chemical potentials. It turns out that in this case there is 
a phase transition at values of /? less than /3c • We have also found a natural 
generalization of the ABC dynamics which leads to a grand canonical Gibbs 
measure: one considers the system on a ring of + 1 sites, with one special 
particle (of type X) replacing the boundaries in the interval model, and 
permits exchanges AX XB, and cyclic permutations of these, at rates 
which satisfy detailed balance with the given chemical potentials. This will 
be the subject of future work. 
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A The profiles as elliptic functions 

The function t{y) defined in (j5.10p is an elliptic function, since Uxiy) is a 
quartic polynomial in y. To identify this function explicitly we introduce a 
fractional linear or Mobius transformation which sends the four roots 0, a, 
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b, and c of Uxiu) to the four roots of the polynomial which appears in the 
incomplete elliptic integral of the first kind [22] , 



F{x,k) 



dt 



y'il-x^Xl-k^x^y 



(A.l) 



We choose the parameter k to satisfy A; > 1 and map (0, a, b, c) 
(— 1, —1/k, 1/k, 1). The formula for the transformation is 



y^z = f{y) 



a+y - 1 

a-y + 1' 



where 



which implies 



±ab + ab{c — 6)(c — a) 



k 



abc 

1 + a_a 



1 — a+a 

Writing the inverse Mobius transformation of ()A.2p as 

1 + z 



z^y = g{z) 



a+ — a_z 



we can write (|5.10p as 



(a+ + a-)dz 



-1/k (a+ - a-zYyJ-2UK{g{z)) 



1/k ^(l-z2)(l-A:2z2) 



where 



2(q!-|_ + a_) 



y^(l — a+c)(l — a+6)(l — a+a) 
The period of oscillation (j5.9p is then 



(A.2) 

(A.3) 
(A.4) 

(A.5) 



(A.6) 



(A.7) 



(A. 
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The inverse of the elhptic integral F{x^ k) is the Jacobi elliptic function 
sn(x, k), that is, F(sn(x, k),k) = x. Then f{y{t)) = s\i{>ct — F /c) , k), so 
that 

a+ — a_ sn(xr — t y-^^k) ,k) 



is the required solution to (jS.Sp . The density profiles are then 

l + sn(2x/32;- |F(i,A:) 



nA{x) 



0+ — a_ sn(2>jr/3x — |F A;) , fc) ' 



nc(a:) 



l + sn(2x/3rg + ^F(i,A;) ,A;) 
0+ — a_ sn(2x/?x + A:) , fc) 



As a check we have derived expressions for values of the elliptic func- 
tions at the special points |F(^,/c) and — |F(^,A;) in a way analogous 
to computing sin57r/12 (i.e., by showing that the desired value is the root 
of a polynomial obtained using trigonometric addition formulae), and then 
used the addition formulae of the elliptic functions [22j to verify directly the 
constancy of the sum (|5.4p and the product (j5.3p of the densities. These 
computations were very intensive and were assisted with the help of com- 
puter algebra package, Maple"""^. 



B Proofs for Remark 15.11 

Proof of Remark \5.1[ (b) It is convenient to change variables in the integrals 
occurring in (j5.9p and (j5.1ip - (j5.12p by writing y = a + {b — a)s, so that for 
example (15. 9p becomes 



TK = 4. . (B.l) 

Jo y/a + s{b — a)y/s{l — s)^yc — {a + s{b — a)) 

Then as K y 1/27 the points a, b, and c approach 1/3, 1/3, and 4/3, 
respectively, so tk approaches 4\/3 ds / \J s(l — s) = 47r\/3- Under this 
change of variable the integral in (jS.lip becomes (assuming that K is so 
small that a < e) 

I ^2 /"^^""^ / ^ + « 

' h ^/s{s + a) V (a + s(& - a))(l - s){c - (a + s{b - a))) ^'^ 
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where e{K) = (e - a)/{b -a). As K \ 0, a = 4K + 0{K^) and b± = 
1 lb 2\J~K + 0(K)^ so that e{K) = e + 0{K) and the second square root m 
()B.2p is r/(s) = 1 + 0{K) uniformly on the integration range, from which 

lim — ^ = hm '^'^ = 2. (B.3) 



K\o ln(l/K) K\o lia{l/K) Jq 7^(7 + 

The argument for the hmit in ()5.12p is similar; the difference in the limiting 
values arises because 1 — 5 is of order ^/K for K small. Finally, a similar 
analysis shows that 

dy 



(B.4) 



\/-2^x(y) 

is bounded as K \ which, with ()B.ip . shows that lim.x\o{TK / ly^i^ / K)) = 
6. 

(c) The proof is long, computational, and unilluminating. We show that 
dTfc/dK > 0; as mentioned in Section \b.2\ it is most convenient to consider 
a G (0,1/3) as the fundamental parameter and show that dTx[a)/da > 
(recall that K{a) = a(l — a)^/4 so dK/da > 0). Our starting point is (|B.ip : 
writing 

/i(a, s) = {a + s{b - a)){c - (a + s{b - a))) (B.5) 



and R = \J a(4 — 3a) we find that 

dTK . f^dh{a,s) ds 4 g{a,s)ds 



da Jo ds /i(a, s)3/Vs(l - ^ Jo h{a,sf/^y/s{l-sy 

(B.6) 

where 

g{a,s) = 3a(l-a)(l+s) + (9a^-12a+2)s2 + (l-3a-(5-9a)s + (2-3a)s^)i?. 

(B.7) 

Unfortunately, g{a, s) is not nonnegative for all a, s, but the sum of g 
and its reflection around the point s = 1/2 is. This follows from the next 
lemma, whose proof is given immediately below. 

Lemma B.l. For < a < 1/3 and < s < 1/2, (i) h{a, s) < h{a, 1 - s); 
(ii) g{a, s) > 0, and (in) g{a, s) + g{a, 1 — s) > 0. 
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Using the lemma, we have from ()B.6P 



da 



> — 



1/2 



+ 



g{a, 1 



g{a,s) 

h{a,sf/'^ ' /i(a,l-s)3/2 
g{a,s) , 5(0,1 



ds 



+ 



0<s<l/2 
g(a,l-s)<0 



/i(a,s)3/2 /i(a,l-s)3/2 



ds 



> y [9{a,s) + g{a,l- s)] 



ds 



(B.8) 



0<s<l/2 
g(a,l-s)<0 



/i(a,s)3/2^s(l -s) 



>0, 



where the three inequahties in (jB.SP are justified by parts (ii), (i), and (iii) 
of the claim, respectively. Note also that the first inequality is strict unless 
g{a, s) < for all s E [1/2, 1], in which case the last inequality will be strict 
(here we use the fact that neither g{a, s) nor g{a, s) + g{a, 1 — s) vanish 
identically in s for any a); thus (jB.SP actually implies drx/da > 0. □ 



Proof of Lemma \B.l[ For (i), note that 

h{a, 1 - s) - h{a, s) = {b - a)(l - 2s){R - a); 



(B.9) 



each of the factors on the right hand side is easily seen to be nonnegative 
for the range of variables in question. Before proceeding to (ii) and (iii) we 
make two preliminary observations. First, the coefficient 



A(a) = (Da^ - 12a + 2) + (2 - 3a)R 



(B.IO) 



of s^ in g{a, s) (which is a quadratic polynomial in s) is positive for all 
a E [0,1/3), for from A(0) = 2, A(l/3) = 0, and A'(l/3) = -8 we see 
that a zero of A in the open interval (0,1/3) would imply two roots of 
A" there, contradicting A'" (a) = 48/i2^ strictly positive. Second, K{a) = 
g{a, 1/2) is positive on [0, 1/3), for K"'{a) = — 6/(a^i?(4 — 3a)) is negative on 
(0, 1/3) and since, as one easily checks, «;"(l/2) = k'(1/2) = k(1/2) = 0, we 
conclude that n"{a), —K'{a), and K{a) are all positive. These observations 
immediately imply (iii) of the claim, since they imply that g{a, s)+g{a, 1 — s) 
is for any a a parabola in s with positive minimum at s = 1/2. 

We turn finally to (ii) of the claim. For fixed a the minimum of g{a, s) 
occurs at 

3a(l - g) + (5 - 9a)R 
- 2(9a2-12a + 2 + (2-3a)i?)- ^^'^^^ 
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For those a for which so(a) > 1/2, the observations in the preceding para- 
graph imply g{a,s) > 0, s G [0, 1/2]. But sq is an increasing function of a, 
smcG 

5{2-3a) + 18a\2-a)+3R 
" ii(9a2 - 12a + 2 + (2 - 3a)i?)2 ' ^^'^^^ 



moreover, if uq = 2/25 then so(ao) = (107V47 - 69)/(88V47 + 686) > 1/2. 
Thus it suffices to show that 

a(-22 + 45a - 36a2 + 9a3) + (4 - 9a + 18a2 - 9a3)i? 

2(9a2-12a + 2 + (2-3a)i?) 

(B.13) 

is nonnegative for a G (0, ao). The denominator in ()B.13p is easily seen to be 
positive in this range. We denote the numerator by (/>(a); then i;^(0) = and 
0(ao) > 0, and positivity will follow once we prove that there is an ai with 
< ai < ao such that (j)'{a) > for a G (0, ai and 4''{a) < for a G (ai, ao). 
But if V(a) = R(l)'{a) then 

tpia) = 108a^-288a3 + 234a^-66a + 8 + (36a^-108a2 + 90a-22)i? (B.14) 

then, as is easily checked, Tp{0) > 0, ip{a) < 0, and ip'{a) < on (0,a + 0), 
since Tp'{a) = {ipi{a) + ip2{0')R)/ where 

V'i(a) = -22 + 303a - SlOa^ + 738a^ - 216a'' (B.15) 

and 

^2(a) = -33 - 432a2 + 234a + 216a^ (B.16) 

when written as polynomials in a — ao, have negative constant term and 
alternating signs, and hence are negative on (0, ao). □ 
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